Theoretical study of Ga-based nanowires and the interaction of Ga with single-wall 

carbon nanotubes 
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Gallium displays physical properties which can make it a potential element to produce metallic 
nanowires and high-conducting interconnects in nanoelectronics. Using first-principles pseudopo- 
tcntial plane method we showed that Ga can form stable metallic linear and zigzag monatomic chain 
structures. The interaction between individual Ga atom and single wall carbon nanotube (SWNT) 
leads to a chemisorption bond involving charge transfer. Doping of SWNT with Ga atom gives 
rise to donor states. Owing to a significant interaction between individual Ga atom and SWNT, 
continuous Ga coverage of the tube can be achieved. Ga nanowires produced by the coating of 
carbon nanotube templates are found to be stable and high-conducting. 

PACS numbers: 73.22.-f, 68.43.Bc, 73.20.Hb, 68.43. Fg, 61.46.+W 



I. INTRODUCTION 

Nanowires display quantum properties which are of in- 
terest from fundamental, as well as technological point 
of view. Quantization of transversal electronic states 
and resulting quantum ballistic conductance and its 
stepwise variation in the course of stretching of wires 
have been studied extensivelyi 1 : 2 : 3 : 4 : 5 : 6 : 7 : 8 : 9 ! 10 The step- 
wise change of the cross section, and l magic l atomic 
structures of a nanowire under tensile stress and the 
interplay between atomic structure and quantized con- 
ductance have been of particular interest in recent 
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From the technological point of view, fabrication of 
stable metallic wires having diameter in the range of 
nanometer has been crucial for the realization of high- 
conducting, low energy dissipating 2 - interconnects in 
nanoelectronics f2i Recent experimental 2 ^ 3 , and theoret- 
ical studies 2 ^ 5 , have demonstrated that such nanowires 
can be produced in a reproducible manner by deposit- 
ing atoms on a single- wall carbon nanotube (SWNT). 
Continuous Ti coating of varying thicknesses, and quasi 
continuous coatings of Ni and Pt were obtained by us- 
ing electron beam evaporation techniques! 22 ' 23 Not only 
metallic connects, but also the contacts of metal elec- 
trodes themselves are crucial for the operation of devices 
based on nanotubesiSi Low resistance Ohmic contacts to 
metallic and semiconducting SWNTs were achieved by Ti 
and Ni atomsiSi The formation of Schottky barrier at the 
contact between metal and SWNT has been found to be 
responsible for the operation of field emission transistors 
made from SWNTs^i 

The electronic and magnetic properties of carbon nan- 
otubes (CNT) can be functionalized by adsorbing dif- 
ferent metal atoms^SiS 9 . It has been shown that transi- 
tion metal atom adsorbed SWNTs have magnetic ground 
states with significant magnetic momenti— The possi- 
bility of filling open nanotubes with liquid by capil- 
lary suction is predicted 3 ^ and filling of the tube with 
molten lead through capillary action is achieved^ More- 



over, temperature measurement by means of a Ga- filled 
CNT thermometer with diameter smaller than 150 nm is 
reported^ Liquid Ga has also potential application as 
micro /nanoswitches. 

Gallium has attracted our interest due to its unusual 
physical properties. It is the only metal, except for mer- 
cury, cesium, and rubidium, which can be liquid near 
room temperatures; this makes its use possible in high- 
temperature thermometers. It has one of the longest liq- 
uid ranges of any metal and has a low vapor pressure 
even at high temperatures. Ga is stable in air and wa- 
ter and high-purity gallium is attacked only slowly by 
mineral acids. Because of these properties Ga can make 
interesting wire structures. Consequently, Ga adsorption 
or Ga coating of SWNT can be of particular interest from 
the electronic devices point of view. 

In this paper we present a systematic study on Ga 
based nanowires. In the first part, binding geometry, co- 
hesive energy, electronic properties and stability of var- 
ious Ga monatomic chain structures are analyzed. In 
the second part, single Ga adsorption on various sites of 
the zigzag (8,0) SWNT is studied. Our prime objective 
is to reveal the character and geometry of the bonding 
and to understand the effect of the adsorption on phys- 
ical properties of SWNT. In this context we also exam- 
ined the formation of a Ga monatomic chain on the (8,0) 
tube. Finally in the third part, we considered the possi- 
bility of producing conducting wires through Ga coating 
of SWNT templates and hence investigated the proper- 
ties of the Ga covered (8,0) tubes. We found that Ga 
monatomic chains are metallic and stable even at 800°K. 
Single Ga atom can be adsorbed to SWNT both from 
outside and inside with a significant binding energy. We 
showed that Ga atoms may form continuous and stable 
coverage of SWNT that transforms the semiconducting 
tube into a high-conducting metal. Our results suggest 
that Ga is an element which can functionalize SWNT to 
make interconnects or contacts for device applications. 
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II. METHOD OF CALCULATIONS 



First-principles total energy and electronic structure 
calculations have been performed using the pseudopo- 
tential plane wave method^within the generalized gra- 
dient approximation (GGA)^ Various structures are 
treated with a periodically repeating tetragonal super- 
cell method. Brillouin-zone integrations are performed 
with 11 special k-points within Monkhorst-Pack special 
k-point schemed All atomic positions in the supercell as 
well as the supercell lattice parameter, c are fully relaxed 
using the conjugate gradient method. The analysis of the 
stability at finite temperature have been carried out by 
relaxing the optimized structures at 800 — 1000°K using 
ab initio molecular dynamics (MD) method with Nose 
thermostat. 

The binding energy (or cohesive energy) Ef, for the 
chain structures is calculated by 



E b = E T [Ga] - E T [Ga - Chain] 



(1) 



in terms of the total energy of single, free Ga atom 
Et [Ga] , and the total energy of chain structure per atom, 
Et [Ga-Chain] , which have been calculated in the same 
supercell. In the study of Ga covered SWNT, the aver- 
age binding energy of Ga atom adsorbed on a SWNT is 
calculated by 

M = {E T [Ga] + E T [SWNT] - Et [SWNT + NGa]}/N 

(2) 

where Et [SWNT] is the total energy (per cell) of fully 
relaxed bare SWNT, and E T [SWNT + NGa] is the total 
energy (per cell) of the SWNT covered by N Ga atoms. 
In these calculations Eb > indicates stable binding. We 
considered the zigzag (8,0) tube as a prototype SWNT 
in our study. The lattice parameter of the bare, relaxed 
(8,0) tube is specified as cswnt = 4.2A, and the cohesive 
energy is calculated to be Eb — 9.1eV per carbon atom. 



III. GA MONATOMIC CHAIN STRUCTURES 

Earlier studies showed that atoms of some 
elements can form stable monatomic chain 
structures! 12 ! 13 ^ 4 ! 18 ! 19 ! 36 ^ 7 ! 38 ^ 9 For example, it has 
been found that Au and Al can form three different 
chain structures j 3 ^ 3 ^ namely the linear chain (LC), 
planar wide angle (WZ) and narrow angle (NZ) zigzag 
structures described by insets in Fig. ^ While LC 
has lowest, NZ has the highest cohesive energy per 
atom. Similar trends have been found also for Group 
IV elements, such as Si, Gej 4 ^ but not C. Interestingly, 
C atom forms only stable LC structure, WZ and NZ 
structures are unstable and change into LC. Usually the 
cohesive energy of LC and WZ are very close. Here we 
analyzed LC, WZ and NZ chain structures of Ga atom. 
Our results are summarized in Fig. ^ The binding 
energy versus lattice parameter, c curves, i.e. Eb(c) per 
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FIG. 1: Variation of binding energy Et of monatomic chain 
structures as a function of lattice parameter, c. Various chain 
structures, linear chain (LC), planar wide angle zigzag (WZ) 
and narrow angle zigzag (NZ) structures are described by in- 
sets. 



atom agrees well with previous first-principles studies 
for Au, Al wires and AuZn, AuMg alloysi 36 i 38 i 39 i 41 
We found among three optimized chain structures NZ 
structure with c = 2.6 A and a ~ 55° has highest 
binding energy (Et, = 2.2 eV), WZ appears a weak 
minimum at c = 4.6 A corresponding to a ~ 131° 
with an intermediate binding energy (Eb = 1.8 eV). LC 
structure has the lowest binding energy (Eb = 1.6 eV) 
with optimized c = 2.6 A and a — 0. 

The stability of these chain structures have been tested 
by raising the temperature to 800°K. To this end we car- 
ried out ab initio MD calculations for 250 time steps 
using a relatively larger supercell containing 4 unit cells. 
We found that LC maintained its string shape, but its 
strict linearity has been destroyed due to the displace- 
ments of atoms at finite temperature. Overall structure 
of NZ also has been maintained, but WZ has been ex- 
cited at high-temperature and eventually trapped at the 
minimum of LC structure. 

Analysis of the electronic band structure shows that 
all three chain structures studied here are metallic. The 
calculated energy band structures are illustrated in the 
Fig. |2 for each geometry. The lowest band of LC is 
due to the er-bond and is derived 4s and 4p z valence 
orbitals of Ga. The doubly degenerate 7r-band is formed 
from the bonding combination of 4p x (4py) orbitals which 
are perpendicular to the axis of LC. Small displacement 
of every second Ga atom perpendicular to the axis of 
the WZ is a distortion that doubles the unit cell but 
lowers the binding energy. The bands of WZ structure 
can be derived from those of LC in Fig. [5J by going from 
LC to WZ all bands of LC are folded and due to the 
breaking of cylindrical symmetry the doubly degenerate 
7T- and 7r*-bands are split. Also the cr*-band of LC which 
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FIG. 2: Left panels: Calculated energy band structures of 
Ga LC, NZ and WZ monatomic chain structures. Zero of the 
energy is set at the Fermi level and shown by dash-dotted 
line. Right panels: Counter plots of the band charge density 
(Pct,Ptt) of LC structure and total charge density (pr) of all 
three chain structures on a plane passing through the Ga-Ga 
bonds. 



dips into the Fermi level near Z, is lifted and becomes 
unoccupied in WZ. The bands of NZ structure are rather 
different from those of LC and WZ structure, since each 
atom of NZ has four nearest neighbors. The character of 
the bonds are illustrated by the charge density contour 
plots in Fig. [21 Except derealization due to cr*-band 
dipping into the Fermi level, the charge density pt{?) 
of LC structure depicts a double bond character. The 
double bond character is strengthen in the WZ structure 
because of the lifting of the <r*-band from the Fermi level. 
The total charge density of the NZ structure becomes 
more uniform and metal-like in the plane of Ga atoms. 

The electrical current through a metallic infinite 
nanowire with infinite mean free path l m — > oo is given 
by Landauer type expression- 2 I = Yjilr\ievi\D(Ep + 
eVh) — Di(Ep)] in terms of degeneracy rji, group velocity 
Vi and density of states of each subband % crossing the 
Fermi levels Then the ballistic conductance becomes 
G = Ei?7i2e 2 /h; namely the number of bands crossing 
the Fermi level times Go (= 2e 2 /h). Based on the band 
structures presented in Fig. [21 the ballistic conductances 
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FIG. 3: A schematic description of different binding sites of 
individual atoms adsorbed on a zigzag (8,0) tube. H: hollow; 
B: bridge; Z: zigzag; T: top site. 



of infinite Ga chain structures are found to be 3Go, 5Go 
and 2Go for LC, NZ and WZ structures, respectively. 
Owing to the small displacement of every second atom 
in WZ structure one channel of LC is closed. In spite 
of two strands in NZ its conductance is smaller than the 
conductance of two parallel LC. 



IV. GA ATOM ADSORPTION ON (8,0) SWNT 

Earlier it has been shown that the electronic proper- 
ties of a SWNT can be dramatically modified by the ad- 
sorption of foreign atomsi 2 ^^. For the adsorption on 
SWNT, we considered four possible sites as an initial po- 
sition of Ga before structure optimization; namely H-site 
above (outside-exo) or below (inside-endo) the surface 
hexagon, the Z- and B-site above the zigzag, and axial 
C-C bonds, and the T-site above carbon atoms. The low- 
est energy binding site is determined by minimizing the 
total energy of SWNT having a Ga atom adsorbed at 
one of those sites described in Fig. [21 To eliminate the 
Ga-Ga interaction, we considered a single Ga atom ad- 
sorbed in every two unit cells of SWNT. Accordingly the 
lattice parameter of the supercell is approximately twice 
the lattice parameter of SWNT, i.e. c ~ 2cswnt- To 
find the lowest binding energy for a given adsorption site 
all atomic positions (SWNT+Ga) and the supercell lat- 
tice parameter have been relaxed. The long range van der 
Waals interaction, Eydw, is expected to be much smaller 
than the chemisorption binding energy and is omitted in 
the present calculations. The calculated E^s for both 
SWNT and graphene are listed in Table |l| The H-site 
is found to be energetically most favorable site among 
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H(exo) 


H(endo) 


B 


T 


Z 


Graphene H 


E b (eV) 


1.7 


2.9 


1.4 


1.4 


1.4 


1.1 


d(A) 


2.4 


2.7 


2.4 


2.3 


2.4 


2.9 



TABLE I: Binding energy Eb and average Ga-C distance of 
Ga adsorbed at different sites on the (8,0) tube. 
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FIG. 4: Energy band structures are calculated for Ga ad- 
sorbed at the H-site of the (8,0) tube with periodicity c ~ 
2cswnt and local density of states, LDOS calculated at 
SWNT and at the adsorbed Ga atom. Upper panels are for 
exo adsorption of Ga; lower ones are for endo adsorption. 
Zero of the energy is set at the Fermi level and shown by 
dash-dotted lines. 



external sites. However, the binding of Ga adsorbed at 
the H-site, but inside the tube yields Eb = 2.9 eV, which 
is stronger than corresponding external adsorption. The 
binding of Ga adsorbed inside is stronger, since the ad- 
sorbate interacts with more C atoms as compared with 
Ga adsorbed outside. According to the Mulliken analy- 
sis 2.1 electrons are transferred from inside Ga to SWNT, 
while only 1.1 electrons are transferred to SWNT when 
Ga is adsorbed outside. These values are consistent with 
Eb l s calculated for inside (endo) and outside (exo) ad- 
sorption. The binding energy is reduced to 1.1 eV for 
Ga adsorbed on the graphene. This is an expected result 
and can be explained by the curvature effect which is 
the primary factor that strengthens the bonding on the 
SWNT. 44 - 4 - 5 - 46 

The electronic structure analysis indicates that Ga ad- 
sorption both inside and outside (in fact, forming a Ga 
LC attached to SWNT with c ~ 2cswnt) metallizes 
the semiconducting SWNT. The calculated energy band 
structures and LDOS for both Ga and SWNT are pre- 
sented in Fig. 01 According to LDOS analysis, the band 
which crosses Ep is derived mainly from SWNT con- 
duction band. This result combined with the Mulliken 






FIG. 5: Counter plots of the total charge density pr(r) of 
exo adsorbed Ga and endo adsorbed Ga on the H-site of the 
SWNT surface. The chemical bonding and charge transfer 
from Ga to SWNT is revealed from counter plots. The axis 
of the SWNT is shown by dashed lines. 



analysis, yielding a charge transfer of 1.1 (2.1) electrons 
from exo (endo) adsorbed Ga to SWNT, suggests that 
Ga electrons are donated to the lowest conduction bands 
of the bare (8,0) tube. As a result, initially empty con- 
duction band of SWNT is gradually populated and also 
modified upon adsorption of Ga, which makes periodic 
Ga+SWNT system (as treated within supercell geome- 
try) a metal. Present results are consistent with the pre- 
vious Al+SWNT studies, 2 ® These results obtained from 
the supercell structure, where Ga atom adsorbed on the 
SWNT is repeated periodically, suggest that the doping 
of a SWNT by an individual Ga atom results in a donor 
state. The charge density contour plots calculated on 
a plane passing through exo (endo) adsorbed Ga atom 
are presented in Fig. [SJ The distribution of charge den- 
sity between Ga and SWNT is consistent with the above 
arguments related with the bonding of Ga atom. 



V. GA ZIGZAG CHAIN ON SWNT 

The discussion in the previous section shows that the 
semiconducting SWNT is metallized as a result of Ga 
adsorption, forming a chain with the lattice parameter 
2cswnt within the supercell geometry. Now we explore 
the effect of the chain formation of adsorbates, and inves- 
tigate the Ga zigzag chain on SWNT. The initial struc- 
ture is obtained by placing a Ga atom at the adjacent 
H-site along the axis the tube. The structure of relaxed 
Ga zigzag chain on SWNT is shown in Fig. HJ1 Here 
the Ga-Ga nearest neighbor distance of 2.6A is close to 
what is obtained in free WZ structure. However, in the 
present case Ga-SWNT distance is increased from 2.4 A 
(the value corresponding to a single exo Ga atom ad- 
sorption) to 2.9A — 3.5A The strong Ga-Ga coupling 
in the adsorbed zigzag chain reduces the charge transfer 
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FIG. 6: Calculated energy band structure, total density of 
states (TDOS) and local density of states (LDOS) at the 
SWNT and at the Ga atoms of the Ga-zigzag chain adsorbed 
on the (8,0) tube. Zero of the energy is set at the Fermi level 
and is shown by dash-dotted line. The adsorption geometry 
is described by inset, where dark and light balls indicate Ga 
and C atoms, respectively. 
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FIG. 7: Intermediate coating obtained by 16 Ga atoms lo- 
cated above each hexagon on the surface of SWNT. Distri- 
bution of Ga-Ga, Ga-C and C-C distances are illustrated by 
gray, dotted and dark histograms respectively in the plot. Ga 
stripes and SWNT template are shown by dark and light balls. 



from Ga atom to the SWNT. This way, while Ga-SWNT 
bond is weakened, the Ga-Ga double bond is formed. 
This situation is consistent with the calculated binding 
energy and charge transfer. Mulliken analysis yields 0.4 
electrons transferred from each Ga atom of the adsorbed 
WZ chain to SWNT. This value is 0.7 electrons smaller 
than that transferred in single Ga adsorption. As for 
the binding energy between zigzag chain and SWNT is 
calculated to be 0.4 eV. 

The strong bonding between Ga-Ga, but weak bonding 
between Ga-SWNT is reflected to the calculated band 
structure shown in Fig. HJ1 Two bands around Ep, one 
is almost fully occupied and the other almost empty are 
reminiscent of the tt- bands of the free WZ in Fig. The 
weak Ga-SWNT interaction causes to distortion of these 
7r-bands. Calculated LDOS's are in compliance with the 
present explanation. 



VI. GA COVERAGE OF SWNT 

Following the adsorption of Ga atom (c = 2cswnt) 
and Ga-zigzag chain (c = cswnt) we next study the Ga 
coverage of SWNT. To this end we attach one Ga atom 
to each H-site on the SWNT; namely 16 Ga atoms per 
unit cell. Upon relaxation of the structure Ga atoms 



move from their original position so that their distance 
from the surface of the SWNT increases. Owing to the 
relatively strong Ga-Ga interaction, Ga-Ga distance is 
reduced while Ga-SWNT distance is increased. At the 
end, Ga islands or stripes are formed above the SWNT 
as shown in Fig. I n the same figure the distribu- 
tion of Ga-Ga, Ga-C and C-C interatomic distances are 
presented by histograms. We note that there are two dif- 
ferent C-C bond distances in the SWNT template. Each 
of these islands or stripes form a metallic domain on the 
surface of SWNT. In order to obtain a more uniform 
coverage we first add four more Ga atoms between four 
islands seen in Fig[7| Upon relaxation Ga atoms rear- 
range and start to form zigzag chain structures on the 
surface of SWNT. 

Finally, two more Ga atoms are added to fill vacant 
sites of Ga coating. Further optimization of the struc- 
ture has led to a continuous coverage with stable zigzag 
Ga structures on SWNT as illustrated in FigBD The 
distribution of interatomic distances is also presented in 
the same figure. In the final structure Ga nearest neigh- 
bor distances range between 2. 5-3.0 A and Ga-C distance 
ranges between 2. 4-4. OA The circular cross section of 
SWNT becomes slightly elliptic. Owing to the weak in- 
teraction between SWNT and Ga, the length of the C-C 
bond and the second nearest neighbor distances do not 
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FIG. 8: Top and side view of fully Ga covered (8,0) SWNT. 
The supercell incorporates 32 C (one unitcell of (8,0) SWNT) 
and 22 Ga atoms. Distributions of first and second nearest 
neighbor Ga-Ga, Ga-C and C-C distances are presented by 
gray, dotted and dark histograms, respectively. Total density 
of states (TDOS), and local density of states (LDOS) calcu- 
lated at the SWNT and at the Ga coating are shown in the 
lower panels. 



display significant changes and variations as compared 
to the bare SWNT. However, the elliptic deformation of 
SWNT gives rise to the dispersions in the third and forth 
nearest neighbor distances. Non-uniformities of the Ga 
coating are attributed to the lattice mismatch between 
SWNT and Ga 2D lattice. The charge transfer is about 
0.4 electrons per Ga atom and like the intermediate cov- 
erage, the system is metallic. It is expected that upon the 
formation of thick Ga coating involving several Ga lay- 
ers around SWNT the charge transfer further decreases 
and Ga-C distance increases. Under these circumstances 
a potential barrier can develop between the metallic Ga 
coating and semiconducting SWNT. Such a situation can 
be exploited as metal-semiconductor junction device. 

The electronic structure having several bands crossing 
the Fermi level indicates that Ga coated SWNT is high- 
conducting. The origin of the metallicity of Ga covered 
SWNT, which was a semiconductor in the absence of 
Ga atoms, is determined by TDOS and LDOS indicating 
that the high state density originate from the Ga atoms. 



The non-uniformity of the interatomic distances which 
may destroy the one-dimensional periodicity can give rise 
to the localization of current transporting states^ This 
is characterized by the localization length £. Owing to 
the tubular nature and infinite cross section we expect 
that for Ga-covered SWNT £ - l m {D/X F ) a with 1< a < 
2. Here D is the diameter of the wire and Xf is the Fermi 
wavelength. Our estimation suggests that £ is in fact 
much larger than typical length scale of interconnects. 

The stability of Ga coating is important to produce 
a stable conducting nanowire. The ab initio MD calcu- 
lation of fully Ga covered SWNT at 1200K has shown 
that the system were stable after 50 time steps. Con- 
tinuity were maintained; neither clusterings nor vacan- 
cies have occurred. Furthermore, we tested the stabil- 
ity of Ga coating (tubular structure) by deleting SWNT 
and by relaxing it at T=0. After full relaxation tubular 
form has remained. It appears that Ga tube made of 
zigzag chains contributes to the stability of Ga covered 
SWNT. The wetting of liquid metals on solid metals and 
non-metallic materials had been studied extensively to 
form composites with exceptional properties^ We be- 
lieve that present results are relevant to the wettability of 
graphite and SWNT by liquid Ga to produce Ga-matrix 
composites or to form Ga intermediates to attach other 
metals. 



VII. CONCLUSION 

The present study has revealed interesting features of 
Ga element in forming chain structure, and also in coat- 
ing of a semiconducting SWNT. Ga forms monatomic 
linear and zigzag chain structures which are metallic and 
stable even at high temperatures. The zigzag structures 
are energetically more favorable than the linear structure. 

The electronic properties of SWNT can be modified 
via single Ga adsorption. Individual Ga adsorbed inside 
or outside the tube gives rise to donor states. This, in 
turn, makes the initially semiconducting SWNT metal- 
lic at finite temperature due to charge transfer from Ga 
atoms to nanotube. SWNT can serve as a template for 
constructing stable Ga chains on its surface. The 4p elec- 
trons from Ga chain partially transferred to SWNT and 
generate additional states around Ep. The band struc- 
ture analysis shows that initially semiconducting SWNT 
transforms into a metal with adsorption of Ga atoms. 

Finally, it is shown that the formation of continuous 
Ga coverage of SWNT can be achieved. This possibil- 
ity revealed by a theoretical study requires, of-course, 
the elaboration of growth conditions which determines 
the quality of growth. The circular cross-section changes 
into elliptical form after the full coverage of Ga. The 
SWNT becomes a quasi-lD metal with high state den- 
sity at Ep. As 3d orbitals of Ga behave like semi-core 
states, the ground state of Ga covered SWNT system is 
non-magnetic. The ab initio MD calculations at 1200°K 
have shown that Ga coating is stable even at high tern- 
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peratures. These results suggest that a SWNT can be 
used as template to produce high-conducting nanowires 
by using Ga coverage. Ga atoms being liquid at room 
temperature and stable in air and water can be used 
to produce ID conductors, thin metallic connects and 
nanometersized electronic devices which may find tech- 
nological applications in nanoscience. In particular, the 
conductivity of liquid Ga-filled SWNT can be exploited 
as nanoswitch. 
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